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A search is presented for the decay T(IS') -7 7 A 0 , A 0 —> cc, where A° is a candidate for the CP-odd 
Higgs boson of the next-to-minimal supersymmetric standard model. The search is based on data 
collected with the BABAR detector at the T(2S) resonance. A sample of T(IS') mesons is selected 
via the decay T(2S) -A n + n~T(lS). The A 0 -A cc decay is identified through the reconstruction of 
hadronic D°. D + , and D* ** (2010) + meson decays. No significant signal is observed. The measured 
90% confidence-level upper limits on the product branching fraction B(T(1S) -A 7 A 0 ) x B(A° —> cc) 
range from 7.4 x 10 -5 to 2.4 x 10 ~ 3 for A° masses from 4.00 to 8.95 GeV/c 2 and 9.10 to 9.25 
GeV/c 2 , where the region between 8.95 and 9.10 GeV/c 2 is excluded because of background from 
T(2S) -7 7Xw(LP), XW(1P) -»• 7^(15) decays. 

PACS numbers: 12.15.Ji, 12.60.Fr, 13.20.Gd, 14.80.Da 


The next-to-minimal supersymmetric standard model 
(NMSSM) is an appealing extension of the standard 
model (SM). It solves the /i-problem of the minimal su¬ 
persymmetric standard model and the hierarchy problem 
of the SM lj, [Ij. The NMSSM has a rich Higgs sector 
of two charged, three neutral CP-even, and two neutral 
CP-odd bosons. Although the Higgs boson discovered at 
the CERN LHC {1, @| is consistent with the SM Higgs 
boson, it can also be interpreted as one of the heavier 
Higgs bosons of the NMSSM Q- The least heavy of the 
NMSSM Higgs bosons, denoted A 0 , could be light enough 
to be produced in the decay of an T meson [1, 6 ). 

In the context of type I or type II two-Higgs-doublet 
models, the branching fractions of the A° depend on the 
A 0 mass and the NMSSM parameter tan/3 Q- Below 
the charm mass threshold, the A 0 preferentially decays 
into two gluons if tan/3 is of order 1 , and to ss or to 
H + fi~ if tan/3 is of order 10. Above the charm mass 
threshold, the A 0 decays mainly to cc for tan /3 of order 
1 and to t + t~ for tan/3 of order 10. BABAR has already 
ruled out much of the NMSSM parameter space for A 0 
masses below the charm mass threshold through searches 
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for and for A 0 —> gg or ss 0 . Above 

the charm mass threshold, BABAR has ruled out some 
of the parameter space for high tan/3 with the A 0 —»• 
t + t~ searches |llUl2|- None of the searches from BABAR 
have observed a significant signal, nor have the searches 
in leptonic channels from the CMS and CLEO EMI 
Collaborations. The A 0 — > cc channel is one of the last 
channels that has not yet been explored. 

We report a search for the decay T(1S) -7 7 A 0 , 
A 0 -7 cc with A 0 masses ranging between 4.00 and 
9.25 GeV/c 2 . An T(1S) decay is tagged by the presence 
of a pion pair from T(2S) -7 An A 0 —> cc 

decay is tagged by the presence of at least one charmed 
meson such as a D°, a D + , or a £>*(2010) + . There¬ 
fore candidates are constructed from the combination of 
a photon, a D meson, and a dipion candidate. An exclu¬ 
sive reconstruction of the A 0 is not attempted. Instead, a 
search is performed in the spectrum of the invariant mass 
of the system that recoils against the dipion-photon sys¬ 
tem. The analysis is therefore sensitive to the production 
of any charm resonance produced in the radiative decays 
of the T(1S) meson. 

The data were recorded with the BABAR detector at the 
PEP-II asymmetric-energy e + e“ collider at the SLAC 
National Accelerator Laboratory. The BABAR detector is 
described in detail elsewhere 0 , 113. We use 13.6 fb " 1 
of “on-resonance” data collected at the T(2S ) resonance, 
corresponding to (98.3 ± 0.9) x 10 6 Y(2S) mesons [l 8 | . 
which includes an estimated (17.5 ± 0.3) x 10 6 T(2S) -A 
tt + tt~T(1S) decays 0. The non-X(25') backgrounds are 
studied using 1.4 fb -1 of “off-resonance” data collected 
30 MeV below the T(2S) resonance. 

The EvtGen event generator 0 is used to simu¬ 
late the signal event decay chain, e + e _ -7- T(2S) -A 
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n +y K~Y(lS), T(1S) —> 7 A 0 , A 0 cc, for A 0 masses be¬ 
tween 4.0 and 9.0 GeV/c 2 in 0.5 GeV/c 2 steps and for A 0 
masses of 9.2, 9.3, and 9.4 GeV/c 2 . The A 0 decay width 
is assumed to be 1 MeV. The hadronization of the cc 
system is simulated using the Jetset [21] program. The 


22 ] 


K®ir + ir , 


detector response is simulated with the GEANT4 
suite of programs. 

Photon candidates are required to have an energy 
greater than 30 MeV and a Zernike moment A 42 [23 
less than 0.1. The A 42 selection reduces contributions 
from hadronic showers identified as photons. Events are 
required to contain at least one photon candidate. Each 
photon candidate is taken in turn to represent the ra¬ 
diative photon in the T(1S) —> 7 A 0 decays. We do not 
select a best signal candidate, neither for the radiative 
photon nor for the D meson and dipion candidates dis¬ 
cussed below, but rather allow multiple candidates in an 
event. 

Events must contain at least one D meson candidate, 
which is reconstructed in five channels: D° —> K~n + 1 
D + —> K~ir + n + , D° —> K~n + n + TT ~, D 1 
and D*(2010) + —> n + D° with D° 

The D° —> K~tt + tt° decays are reconstructed in the 
D*( 2010 ) + production channel to reduce a large back¬ 
ground that would otherwise be present. The inclusion 
of charge conjugate processes is implied. The 7 r° candi¬ 
dates are reconstructed from two photon candidates by 
requiring the invariant mass of the reconstructed 7r° to 
lie between 100 and 160 MeV/c 2 . The 7 r° candidates do 
not make use of the radiative photon candidate. The I\ ° 
candidates are reconstructed from two oppositely charged 
pion candidates. Each I\® candidate must have a re¬ 
constructed mass within 25 MeV/c 2 of the nominal K® 
mass EH and satisfy d/ad > 3, where d is the distance 
between the reconstructed e + e _ collision point and the 
K® vertex, with ad the uncertainty of d. 

The D° and D + candidates are required to have masses 
within 20 MeV/c 2 of their nominal masses (Hi corre¬ 
sponding to 3 to 4 standard deviations (a) in their mass 
resolution. When reconstructing H*(2010) + candidates, 
we constrain the D° —> K~tt + tt 0 candidate mass to its 
nominal value EH- The Z?*(2010) + candidate mass dis¬ 
tribution has longer tails. The D*(2010) + candidates are 
required to lie within 5 MeV/c 2 of its nominal mass EH, 
corresponding to 10 a in the mass resolution. 

Events are required to have at least one dipion can¬ 
didate, constructed from two oppositely charged tracks. 
The invariant mass, m Rl of the system recoiling against 
the dipion in the Y[2S) —> tt + tt~Y(1S) transition is cal¬ 
culated by 


after application of these criteria. A clear peak is seen at 
the T(1S) mass. 
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FIG. 1: The rriR distribution of events with a dipion, charm, 
and photon tag before application of selection criteria based 
on the BDT output (see text). The solid circles indicate 
the on-resonance data. The open squares indicate the off- 
resonance data normalized to the on-resonance luminosity. 

All charged tracks and calorimeter clusters other than 
those used to define the radiative photon, the D meson 
candidate, and the dipion candidate are referred to as 
the “rest of the event”. 

The mass of the A 0 candidate, mx , is determined from 
the mass of the system recoiling against the dipion and 
photon through 


m 2 X = (P e +e~ ~ Pn+ir- ~ Pjf, 


( 2 ) 


m 2 R — M'y^2S) + m T7i- — 2M r (2S).E 7r7r , 


(1) 


where P denotes four-momentum measured in the CM 
frame. The four-momentum of the e + e _ system is given 
by Pe+e- = {Mt(2S), 0,0,0). 

Backgrounds are evaluated using simulated T(2S) and 
e + e _ —► qq events, where q is a u, d , s, or c quark. Events 
with low-energy photons contribute a large background 
for mx greater than 7.50 GeV/c 2 . Therefore, the analysis 
is divided into a low A 0 mass region (4.00 to 8.00 GeV/c 2 ) 
and a high A 0 mass region (7.50 to 9.25 GeV/c 2 ). The 
definitions of the regions, which overlap, are motivated 
by the need to have sufficient statistical precision for the 
background determination in each region. 

We train ten boosted decision tree (BDT) classi¬ 
fiers [ 3 ] to separate background from signal candidates 
(two mass regions x five D channels). The BDTs are 
trained using samples of simulated signal events, simu¬ 
lated generic Y(2S) events, and the off-resonance data. 
The BDT inputs consist of 24 variables: 

1-2. Event variables: 


where m W7r is the measured dipion mass, My^S) is the 
nominal T(2S') mass [l9j|> an d P-mr is the dipion energy in 
the e + e _ center-of-mass (CM) frame. The two pions in 
the dipion system are required to arise from a common 
vertex. Signal candidates must satisfy 9.45 < m R < 
9.47 GeV/c 2 . Figure [T] presents the distribution of mu 


• number of charged tracks in the event, 

• number of calorimeter clusters in the event. 

3-12. Kinematic variables: 

• rn R , 
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• dipion likelihood (defined later), 

• D candidate mass, 

• D candidate momentum, 

• photon 7T° score (defined later), 

• energy of the most energetic charged track in 
the rest of the event, calculated using a charged 
pion mass hypothesis, 

• energy of the most energetic calorimeter cluster 
in the rest of the event, 

• invariant mass of the rest of the event, 

• CM frame momentum of the rest of the event, 

• CM frame energy of the rest of the event. 

13-15. Vertex variables: 

• transverse coordinate of a vertex formed using 
all charged tracks, 

• longitudinal coordinate of a vertex formed using 
all charged tracks, 

• the x 2 probability of a vertex fit using all 
charged tracks. 

16-18. Event shape variables: 

• the ratio of the second to zeroth Fox-Wolfram 
moment HH , calculated using all charged tracks 
and calorimeter clusters, 

• sphericity [ 2 f| of the event, 

• magnitude of the thrust [ 2 T| . 

19-24. Opening angles in the CM frame between the 

• dipion and photon candidate, 

• dipion and D candidate, 

• dipion and thrust axis, 

• photon and D candidate, 

• photon and thrust axis, 

• D candidate and thrust axis. 

The kinematic variables provide the most separation 
power for all ten BDTs. The separation power of the 
other variables depends on the mass region and channel. 
The vertex variables suppress background without a D 
meson in the event. The event shape variables suppress 
e + e _ —> qq backgrounds. 

The dipion likelihood 0 is defined using the opening 
angle between the two charged pions in the CM frame, 
the transverse momentum of the pair, the invariant mass 
of the pair, the larger of the two momenta of the pair, 
and the y 2 probability of the pair’s vertex fit. 

To reject photon candidates from tt° —» 77 decays, 
a likelihood [24j is defined using the invariant mass of 
the radiative photon candidate and a second photon (if 
present), and the second photon’s CM energy. The lower 
the likelihood, the more 7 r°-like the photon pair. The 


photon 7T° score is the minimum likelihood formed from 
the radiative photon and any other photon in the event 
excluding photon candidates used to reconstruct the 7r° 
candidate in the D° —> decay. 

For each channel and mass range, each BDT output 
variable is required to exceed a value determined by max¬ 
imizing the quantity 5/(1.5 + \f~B) [28), where 5 and 
B are the expected numbers of signal and background 
events, respectively, based on simulation. 

In the case of events with multiple signal candidates 
that satisfy the selection criteria, there may be multiple 
values of mx ■ Signal candidates that have the same dip¬ 
ion and radiative photon candidate have the same value 
of uix , irrespective of which D candidate is used. We 
reject a signal candidate if its value of mx has already 
been used. 

In total, 9.8 x 10 3 and 7.4 x 10 6 candidates satisfy the 
selection criteria in the low- and high-mass regions, re¬ 
spectively. The corresponding distributions of mx are 
shown in Fig. [2] The backgrounds in the low-mass re¬ 
gion consist of T(15) — > 7 gg (35%); other T(15) de¬ 
cays, denoted T(15) — > X (34%); T(25) decays with¬ 
out a dipion transition, denoted T(25) —> X (15%); and 
e + e _ —> qq events (16%). The corresponding background 
contributions in the high-mass region are 1 %, 66 %, 18%, 
and 15%. Background contributions from T(15) —► 7 gg 
decays reach a maximum near 5.5GeV/c 2 and decrease 
above 7 GeV/c 2 . 

We search for the A 0 resonance as a peak in the mx 
distribution. The reconstructed width of the A 0 is ex¬ 
pected to strongly depend on its mass due to better pho¬ 
ton energy resolution at lower photon energies. There¬ 
fore, an extended maximum likelihood fit in a local mass 
range is performed as a function of test-mass values, de¬ 
noted mA°- For these fits, the parameters of the proba¬ 
bility density function (PDF) used to model the shape of 
the signal distribution are fixed. The parameters of the 
background PDF, the number of signal events N S i g , and 
the number of background events are determined in the 
fit. 

The si gna l mx PDF is modeled with a Crystal Ball 
function [29j) . which consists of a Gaussian and a power- 
law component. The values of the signal PDF at a given 
value of m.A° are determined through interpolation from 
fits of simulated signal events at neighboring masses. The 
background mx PDF is modeled with a second-order 
polynomial. 

The fits are performed to the mx spectrum, for vari¬ 
ous choices of m^o, in steps of 10 and 2 MeV/c 2 for the 
low- and high-mass regions, respectively. The step sizes 
are at least 3 times smaller than the width of the sig¬ 
nal mx PDF. We use a local fitting range of ±10 cr qb 
around m^o, where <jcb denotes the width of the Gaus¬ 
sian component of the Crystal Ball function. The acB 
parameter varies between 120 and 8 MeV/c 2 for values of 
to^o between 4.00 and 9.25 GeV/c 2 , as shown in Fig. [3] 
We do not perform a fit for 8.95 < to^o < 9.10 GeV/c 2 
because of a large background from T[2S) —> 7 Xf,j(lP), 
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FIG. 2: The mx distributions of signal candidates in the 
low- (a) and high- (b) mass regions after applying all selec¬ 
tion criteria. The points indicate the data. The solid curve 
shows the result of a fit to the data under a background- 
only hypothesis. The colored histograms show the cumulative 
background contributions from e + e _ —>• qq (magenta dense- 
dot filled), X(2S) — > X (green sparse-dot filled), T(1S) — > X 
(blue dotted), and T(IS') —» 7 gg (red dashed) events. 


XbjO-P) decays. 

The fitting procedure is validated using background- 
only pseudo-experiments. The rax PDF used to generate 
pseudo-experiments for the low-mass region is obtained 
from a fit of a fifth-order polynomial to the low-mass 
region data. The mx PDF used for the high-mass re¬ 
gion is obtained from a fit of the sum of four exponential 
functions plus six Crystal Ball functions to the high-mass 
region data, with shape parameters fixed according to ex¬ 
pectations from simulation and with the relative normal¬ 
izations determined in the fit. The Crystal Ball functions 
describe the T(2S) -X 7Xbj(lP) and Xbj(l-P) —> S) 

transitions while the exponential terms describe the non¬ 
resonant background. Four exponential terms are used 
because the non-resonant background increases rapidly 
for higher mx- The background fits are overlaid in Fig. [2] 
The fitting procedure returns a null signal for most to ,40 
values but is found to require a correction to N S i g for 
values of near 4.00 or 9.25 GeV/c 2 . The correc- 
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FIG. 3: The <jcb parameter for A 0 decays of various simu¬ 
lated masses. 


5 6 7 8 9 
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tions are determined from the average number of signal 
events found in the fits to the background-only pseudo¬ 
experiments. The corrections are applied as a function 
of to ,40 and reach a maximum of 15 and 50 candidates in 
the low- and high-mass regions, respectively. The uncer¬ 
tainty of the correction is assumed to be half its value. 

The reconstruction efficiency takes into account the 
hadronization of the cc system into D mesons, the 
branching fraction of D mesons to the five decay chan¬ 
nels, detector acceptance, and the BDT selection. The 
efficiencies range from 4.0% to 2.6% for simulated A 0 
masses between 4.00 and 9.25 GeV/c 2 . 

Potential bias introduced by the fitting procedure is 
evaluated using pseudo-experiments with different values 
of the product branching fraction B(T(1S) —> 7 A 0 ) x 
B(A° —> cc). For various choices of 777 , 40 , the extracted 
product branching fraction is found to be (4±3)% higher 
than the value used to generate the events. This result 
is used to define a correction and its uncertainty. 

Table U summarizes all correction factors and associ¬ 
ated systematic uncertainties. The fit correction system¬ 
atic uncertainty is added in quadrature with the statisti¬ 
cal uncertainty of N S i g . All other correction factors are 
added in quadrature and applied to the reconstruction ef¬ 
ficiency. A correction of 1.00 means we do not apply any 
correction but propagate the multiplicative uncertainty. 

The systematic uncertainties associated with the re¬ 
construction efficiencies are dominated by the differences 
between data and simulation, including the BDT out¬ 
put modeling, cc hadronization, D-candidate mass reso¬ 
lution, dipion recoil mass and likelihood modeling, and 
photon reconstruction. Other systematic uncertainties 
are associated with the fit bias (discussed above), the 
dipion branching fraction 0 , the finite size of the sim¬ 
ulated signal sample, and the T(2S) counting fl8j |. 

The BDT output distributions in off-resonance data 
and e + e _ —> qq simulation, shown in Fig. [4] have con¬ 
sistent shapes but are slightly shifted from one another. 
The associated systematic uncertainty is estimated by 
shifting the simulated distributions so that the mean val- 
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TABLE I: Summary of corrections and their associated sys¬ 
tematic uncertainties. All corrections are multiplicative ex¬ 
cept for the fit correction. 


Source 

Low region High region 

Fit correction (candidates) 

up to 15 ± 8 up to 50 ± 25 

BDT output modeling 

0.93 ±0.04 0.98 ±0.01 

Source 

Both regions 

cc hadronization 

1.00 ±0.09 

Fit bias 

1.04 ±0.03 

Dipion branching fraction 

1.00 ±0.02 

Photon reconstruction 

0.967 ±0.017 

D mass resolution 

0.98 ±0.01 

Finite simulation statistics 

1.00 ±0.01 

T(2S) counting 

1.00 ±0.01 

Dipion likelihood 

1.02 ±0.01 

Dipion recoil mass 

0.991 ±0.005 


is estimated by Gaussian smearing of the D mass input 
in simulation to match the data and measuring the dif¬ 
ference in the reconstruction efficiency. 

Further corrections to account for data and simulation 
differences in reconstruction efficiencies are estimated 
with similar methods. Corrections are applied to account 
for the dipion recoil mass reconstruction, the dipion like¬ 
lihood modeling, and the photon reconstruction [iioj . 

The highest observed local significance in the low-mass 
region is 2.3 standard deviations, including statistical un¬ 
certainties only, at 4.145 GeV/c 2 . The corresponding re¬ 
sult for the high-mass region is 2.0 standard deviations 
at 8.411 GeV/c 2 . The fits are shown in Fig. [5] Such fluc¬ 
tuations occur in 54% and 80% of pseudo-experiments, 
respectively. Hence our data are consistent with the 
background-only hypothesis. 


ues agree with the data, and then recalculating the effi¬ 
ciencies. The reconstruction efficiencies decrease by 7% 
and 2% in the low- and high-mass regions, respectively. 



BDT output 


FIG. 4: The BDT distributions in off-resonance data (points) 
and simulated e + e“ — ¥ qq events (histograms) for the five D 
meson decay modes. The results on the left (a, c, e, g, i) 
and right (b, d, f, h, j) correspond to the low- and high-mass 
regions, respectively. 



The uncertainty associated with cc hadronization is 
evaluated by comparing D meson production in off- 
resonance data and e + e _ —> cc simulation normalized 
to the same luminosity. The difference in the yield varies 
from 1% to 9% for the five D decay channels. We conser¬ 
vatively assign a global multiplicative uncertainty of 9% 
that includes effects due to the hadronization modeling, 
particle identification, tracking, n° reconstruction, and 
luminosity determination of the off-resonance data. 

The uncertainty due to the discrepancy between the 
reconstructed D mass resolution in data and simulation 


FIG. 5: The fits with the highest local significance in the 
low- (a) and high- (b) mass regions. The solid line is the 
fit that includes a signal. The dotted line is the background 
component of the solid line. 

Upper limits on the product branching fraction 
B(T(1S) —l 7^4°) x B(A° —> cc) at 90% confidence level 
(C.L.) are determined assuming a uniform prior, with the 
constraint that the product branching fraction be greater 
than zero. The distribution of the likelihood function for 
N S i g is assumed to be Gaussian with a width equal to 









































the total uncertainty in N S i g . The upper limits obtained 
from the low-mass region are combined with those from 
the high-mass region to define a continuous spectrum for 
the upper limits. The results are shown in Fig. [6] 


ing from 7.4 x 10~ 5 to 2.4 x 10 -3 for A 0 masses from 
4.00 to 9.25 GeV/c 2 , excluding masses from 8.95 to 9.10 
GeV/c 2 because of background from T(2S) —> 7Xbj(lP), 
Xbj(lP) —> 7T(15) decays. These results will further 
constrain the NMSSM parameter space. 



FIG. 6: (color online) The 90% C.L. upper limits on the 
product branching fraction B(T(1S) —)■ 7A 0 ) x B(A° —> cc) 
using all uncertainties (thick line) and using statistical uncer¬ 
tainties only (thin dashed line). The inner and outer bands 
contain 68% and 95% of our expected upper limits. The bands 
are calculated using all uncertainties. The thin solid line in 
the center of the inner band is the expected upper limit. 


In summary, we search for a resonance in radiative de¬ 
cays of the T(1S) with a charm tag. We do not observe 
a significant signal and set upper limits on the product 
branching fraction B(T(1S) —> 7 A 0 ) x B(A° —> cc) rang¬ 
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